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A connection between radio state changing and glitch 
activity in PSR J0742-2822 
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ABSTRACT 

PSR J0742— 2822 exhibits two distinct emission states that are identified by discrete 
changes in the observed pulse profile. These changes have previously been shown to 
correlate with changes in the derivative of the pulse frequency. In this paper, we 
use observations with the Parkes radio telescope at a centre frequency of 1369 and 
3100 MHz to produce high phase resolution polarisation profiles for the two modes 
and perform a detailed study of the correlation between observed pulse shape and 
spin-down rate. We find no correlation for at least 200 days prior to a glitch in the 
pulsar at MJD 55022, following which the correlation becomes strong. This suggests a 
link between the emission state switching phenomenon and glitch events. We discuss 
the possibility that emission state switching is driven by the interaction between the 
magnetosphere and the interior of the neutron star. 
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1 INTRODUCTION 

Radio pulsars are renowned for highly regular pulsations 
driven by their stable rotation. Although some pulsars 
have precision of tens of nanoseconds over many years 
IjManchester et alj 12013 ). this is not the case for the ma- 
jority of pulsars, especially those with characteristic ages 
less than ~ 100 kyr. 

There are t wo main d r ivers of instability in pulsars, 
"glitches" (e.g. ILvne et all l2000l: lAtoar et all Il98i) and 
"timi ng noise" (e.g. ICordes fc Helfandl Il98d : iHobbs et al.l 



|2010j). Glitches are characterised by a sudden change in the 
observed pulse period. These sudden changes are typically 
followed by a slower "recovery" , where the period returns to- 
wards the value it would have had without the glitch. This 
recovery occurs over a wide range of time-scales, although 
it is a t present not well understood (|van Evsden fc Melatosl 
l2010l h 

Timing noise describes the unmodelled red noise in 
the observed pulse arrival times of t he majority of ra- 
dio pulsars dCordes fc Downs! 1 19851 : IHobbs et all l20ld : 
IShannon fc Corded |2010j). The physical processes behind 
both of these observationally defined phenomena are not 
well understood, however there are two main classes of the- 
ory. Either the variation in pulse arrival times is driven by 
the interior of the neutron star, or due to a change in the 
magnetospheric arrangement. 

Most accepted theories for glitches model their origin 
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in the interior of the neutron star, where angular momen- 
tum is stored and then suddenly transferred to the crust 
in order to cause a sudden change in the spin of the neu - 
tron star surface (|Anderson fc Itohlll975l : lAlpar et al.lll984T ). 
Timing noi s e is le ss well understood; however recent work by 
ILvne et al.l (|2010j) demonstrated that for six pulsars, the ob- 
served timing noise (characterised by fluctuations in period 
derivative) correlates with the observed pulse shape (charac- 
terised by pulse width or relative amplitude of pulse profile 
components). This suggests that some timing noise may be 
driven by some magnetospheric process, possibly with dis- 
crete state switching on timescales from 100 days to greater 
than 2000 days. Similar profile changes have also been ob- 
served in PSR J0738— 404 2, with the profile shape varying 
on timescales of decades (jKarastergiou et al.ll201lh . It has 
been suggested that the timing noise in the youngest pul- 
sars is dominated by unmodelled reco very from gl itches and 
there fore by the neutron star interior (|Johnston fc Galloway! 
Il999j), but as the pulsar ages, glitch activity decreases and 
the timing noise is dominated by changes in the magneto- 
sphere. However, there are some indications that glitches 
can also be linked to pulse shape changes. A recent glitch 
in PSR J1119— 6227 coincided with the appearance of ad- 
ditional pu lse components with in termittent or RRAT-like 
behaviour jWeltevrede et aLlfeoilT ). 

In this paper, we present a study of PSR J0 742— 2822, 
one of the six pulsars studied bv lLvne et alj 1)20101) . This pul- 
sar exhibits the most rapid state changes of the sample, and 
has the least clear correlation between pulse shape changes 
and spin-down changes. In addition, this pulsar has several 
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published glitches , with the most recent (MJD 55022) being 
of moderate size l|Espinoza et al.ll201ll ; lYu et al.ll2012l ). 



2 OBSERVATIONS 

Since 2007, PSR J0742-2822 has been observed on a 
roughly monthly basis at 1369 MHz and six-monthly at 
3100 MHz with the Par kes radio telescope as p art of the 
Fermi timing programme (|Weltevrede et al.ll201Cl) . The data 
were recorded with one of the Parkes Digital Filterbank sys- 
tems (PDFB1/2/3/4), with a total bandwidth of 256 MHz 
at 1369 MHz and 1024 MHz at 3100 MHz. We supplemented 
these observa tions with archival data from the Parkes pulsar 
data archive (|Hobbs et al.ll201ll ). 

Flux densities have been calibrated by comparison to 
the continuum radio source 3C 218. Observations are cali- 
brated for differential gain and phase between the feeds us- 
ing measurements of a noise diode coupled to the receptors 
in the feed. To correct for cross-coupling of the receptors in 
the Multibeam receiver used for the 1369 MHz observations, 
we used a model of the Jones matrix for the receiver com- 
puted by observation of the bright pulsar PSR J0437— 4715 
over the entire range of hour angles visi ble, using the 'mea - 
surement equation modelling' technique (|van Stratenll2004T) . 
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3 TIMING 

We follo w the same general pu lsar timing procedure as 
used by IWeltevrede et al.l (|2010r ) to obtain a best-fit pul- 
sar ephemeris. Because of the large amount of timing noise, 
we fit t he timing model using the pre-whitening tech- 
nique of I Coles et al.l (|201l[ ) to perform an unbiased least- 
squares fit in the presence of red noise. This model includes 
pulse frequency, frequency derivative, position, proper mo- 
tion and the previously known glitch event at MJD 55022 
l|Espinoza et all 1201 ll ; lYu et al.l I2012T ). however it does not 
attempt to model the timing noise. We measure the step 
change in pulse frequency associated with the glitch to be 
Af = 0.61 fjHz , cons i stent with that previously reported 
l|Espinoza et all 1201 ll ; lYu et alll2012J ). The difference be- 
tween our best-fit model and the observed time of arrival 
is termed the timing residual, and in the case of PSR 
J0742— 2822 is dominated by a combination of red timing 
noise and glitch recovery. 

The pulsar observations are irregularly sampled and so 
we fit a smooth curve to the timin g resid uals using the inter- 
polation technique of iDeng et all (|2012j), which gives us an 
effective time resolution of 60 days. Both the residuals and 
the smooth fit are shown in Figure [T] When the slope of the 
residuals is positive the pulse frequency is overestimated. By 
taking the gradient of the residuals, we can determine the 
anomalous pulse frequency, i.e. the fluctuations which are 
not modelled by our ephemeris. Similarly, the second deriva- 
tive gives us the anomalous frequency derivative. These are 
also shown in Figure [T] 

We confirm that the frequency derivative shows quasi- 
perio dic oscillations w ith a timescale on the order 100 
days (|Lvne et all 120101 ). We also note that the mean fre- 
quency derivative after the glitch is reduced by about 
which we attribute to the step change in 



Figure 1. Timing noise in PSR J0742— 2822. Upper panel shows 
the post-fit timing residual for PSR J0742-2822. Dots show indi- 
vidual observations (with error bars too small to show), to which 
we fit an optimally constrained smooth curve (solid line). These 
residuals have been corrected for the discrete change in pulse fre- 
quency due to a glitch at MJD 55022, marked with a vertical 
dashed line. The middle panel shows the unmodelled variation in 
pulse frequency, derived from the derivative of the smooth fit to 
the residuals. The lower panel shows the derivative of the above, 
the unmodelled variation in pulse frequency derivative. 



v caused by the g litch. Thi s value of Av is consistent with 
that reported by lYu et all (120121) , rath er than the larger 
value reported by lEspinoza et all (|201ll ) . This step in Az> 
has not been removed from the Figure [T] but we correct for 
it in the remainder of this work. 



4 PULSE PROFILES 

Average pulse profiles of PSR J0742-2822 exhibit com- 
plex stru cture wi t h as many as seven overlapping com- 
ponents l|Kramerl 1 19941 ). The pulses are highly linearly 
polarised, with a small amount of negative circular and 
slight depolarisation towards the trailing edge of the 
pulse (|Karastergiou fc Johnstonll20061 ). The viewing geome- 
try is not well constrained, but it is likely to be close to or- 
thogonal and the pr ofile is symmetrically distributed around 
the magnetic pole (I Johnston et al.ll2005l ). The evolution of 
the profile with frequency is constant with a central 'core' 
component surrounded by at least two 'conaP rings, though 
it has been noted that the relative spectral index of some 
components appears to change significantly with frequency 
l| Johnston et alJl2006t ). 

Figure [2] shows a large number of integrated observa- 
tions of the pulsar at a centre frequency of 1369 MHz, taken 
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Figure 2. Stacked profiles of 380 observations of PSR 
J0742-2822 at a centre frequency of 1369 MHz, normalised to 
keep the pulse width constant. Note that the observations are 
not uniformly spaced in time. 



over five years, stacked sequentially. These profiles are aver- 
ages of typically three minutes of observation and have been 
normalised by the leading and trailing edges to account for 
variation in overall flux density. Two emission modes are 
clearly identifiable, with the most significant change occur- 
ring between phase ~ —4° to —1° on this figure. We define 
the most frequently observed mode to be "Mode I" and the 
less frequently observed mode to be "Mode II". The two 
modes are not completely distinct and some profiles have an 
intermediate shape showing characteristics of both modes. 
Nevertheless, we can split the profiles into the two modes 
and compare the high-resolution average profiles. To clas- 
sify the profiles we define a shape parameter based on the 
1369-MHz pulse profile to be the ratio of the first and second 
components in the profile (the two leading peaks in Mode I; 
see Figure [3}. For Mode I, this parameter is roughly unity, 
and for Mode II this is less than unity. We find that this 
method is supported by subjective assessment of individual 
profiles, and although other parameters can be derived we 
find that the choice of parameter does not affect the subse- 
quent analysis. 

The upper two panels of Figure [3] show the average 
profiles for the two modes of PSR J0742-2822 at a cen- 
tre frequency of 1369 MHz. Here we have averaged a total 
of ~ 15000 s for Mode I and ~ 8000 s for Mode II which ren- 
ders day-to-day flux density variations caused by interstellar 
scintillation negligible, and so normalisation is not required. 
The two modes have a very similar overall pulse width, and 
are very similar on both the leading and trailing edge of the 
profile. Mode II exhibits significantly more emission in the 
central part of the profile, with prominent increases in the 
second and third peaks of the profile as well as a general 
increase of the trailing edge. Both modes retain a high lin- 
ear polarisation fraction, however Mode II is slightly more 
polarised on the trailing edge. The position angles are also 
consistent within their measurement uncertainties, except 
for an orthogonal mode jump observed in Mode I at the 
very trailing edge of the profile. 

An identical analysis was performed for the observa- 
tions centred at 3100 MHz, and the results are shown in 
the lower two panels of Figure [3] The change between the 
two modes is even more pronounced at the higher frequency, 
with the Mode I profile dominated by a leading narrow com- 
ponent. The Mode II profile is almost rectangular, with the 
trailing half of the profile almost doubling in flux density 
over Mode I. The Mode I profile again shows an orthogonal 
polarisation mode jump at the trailing edge of the profile, 
though over a different range of pulse phase to that seen 
at 1369 MHz. Although not surprising, we note that the 
Mode I and II profiles are consistent between the two fre- 
quencies, i.e. when the pulsar is in a particular mode at one 
frequency, it is in that same mode at the other frequency. 
The leading component, which dominates Mode I, does not 
appear to change amplitude significantly between the modes 
and has a flatter spectral index than the components that 
dominate in Mode II. We note that the increased ampli- 
tu de of the pulse-sh ape changes at high frequency may have 
led lJohnston et al.l (J2006) to confuse emission mode changes 
with frequency evolution in their multi-frequency study. 

A number of observations were recorded with sufficient 
time resolution to study individual pulses. Figure [4] shows 
the phase resolved modulation index (ratio of variance to 
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Figure 3. Average pulse profiles of PSR J0742-2822 for Mode 
I (black) and Mode II (grey). The upper two panels, (a) and (b), 
are for a centre frequency of 1369 MHz, and the lower two panels, 
(c) and (d) are for a centre frequency of 3100 MHz. Panels (a) 
and (c) show the polarisation position angle as a function of pulse 
phase. In panels (b) and (d), total intensity is shown with a thick 
line, linear polarisation with a thin line and circular polarisation 
with a dotted line. Note that for much of the profile the pulse 
is almost 100% linearly polarised and so the line showing linear 
polarised intensity is obscured by the total intensity. 
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Figure 4. Two observations of PSR J0742-2822 recording 1000 
individual pulses, taken on 2012-04-15 (Mode I; black lines) and 
2012-05-27 (Mode II; grey lines). The solid lines show the modu- 
lation index and dotted lines show the average profile from that 
observation for reference. 



mean of individual pulses) for representative observations in 
each of the two modes. Due to a hardware limitation when 
recording individual pulses, the phase resolution of this data 
is half that in Figure [3] Even though the observations show 
changes in the mean profile corresponding to the two modes, 
the underlying statistics of the individual pulses remain 
cons tant and in agreement with previously published val- 
ues IjWeltevrede et al.ll2006l ). Both leading components have 
similar modulation indices, even though one of these com- 
ponents varies considerably between modes and the other is 
almost unchanged. The central peak seen in the modulation 
index is associated with a narrow component that is seen to 
vary considerably between observations but does not appear 
to be strongly correlated with the overall mode changing or 
timing analysis. 



5 ANALYSIS AND DISCUSSION 

In Figure [5] we overlay the pulse shape parameter defined 
in Section [3] computed for each observation, with the ob- 
served frequency derivative deviations described in Section 
[3] We also show the pulse shape parameter averaged un- 
der a running 60-day window to match the time resolution 
of the ii measurements. The correlation between spin-down 
rate and pulse shape parameter is very clear for much of 
the timespan. We find that our Mode II is associated with 
increased \i>\, i.e. larger negative v. This association of the 
mode dominated by steep spectral index components with a 
larger spin-down rate supports the hypothesis that so-called 
central 'core' co mponents tend to dominate the large spin- 
down rate mode (lLvne et al.l 2010). 

For two years prior to the MJD 55022 glitch of PSR 
J0742— 2822 the pulse profile exhibits two distinct emission 
states, however this does not seem to correlate strongly with 
the pulsar spin down parameter. Figure |S] shows the corre- 
lation between the pulse shape parameter and Av for over- 
lapping 300-day windows. The correlation is computed af- 
ter averaging the shape parameter measurements under a 
running box-car of width 60 days to match the effective 
time resolution of the Ai> measurements. This confirms that 
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Figure 5. Fluctuations in frequency derivative of PSR J0742— 2822 in excess of the best-fit timing model (See Section [3}. Circles 
indicate the measured pulse shape parameter for each of the observations. A grey dotted line shows the shape parameters averaged under 
a running box of width 60 days. The vertical dashed line at MJD 55022 indicates the epoch of a glitch. 
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Figure 6. Correlation between frequency derivative and 
smoothed pulse shape parameter for overlapping 300-day inter- 
vals. The vertical dashed line at MJD 55022 indicates the epoch 
of a glitch. Also shown with dotted bars is the same correlation 
when computed for the entire pre and post-glitch epochs. 



the correlation swiftly increases at the glitch event of MJD 
55022 and remains high for more than 1000 days. Although 
we observe a change in correlation, it can also be possible 
that the rate of state switching is too rapid prior to the 
glitc h for us to detect with our 60-day time resolution. In- 
deed [Lyne_et_alJ (|2010T ) show that the mode-switching rate 
of PSR J0742-2822 is the most rapid of the sample and 
that this rate is not constant over the entire data-span. As 
we only have data for one glitch occurrence this change in 
correlation co-efficient may be coincidental, but we feel that 
it is worthwhile to consider the possibility that the change 
in the state- switching behaviour is linked to the glitch event. 
Glitch events are generally thought to be driven by 
the neutron star interior, but the emission and spin-down 
is thought to be driven by the magnetosphere. There are 
suggestions that a link between pulse shapes and glitches 
is present in PSR J1119— 6227, where an unusual glitch ap- 
peared to trigger the appearance of additional components 



in the profile IjWeltevrede et al.ll201ll ). However, unlike PSR 
J1119— 6227, we do not see any evidence that the change in 
emission in PSR J0742-2822 is associated with RRAT-like 
statistics of individual pulses. We still do not have any clear 
picture of how the magnetospheric changes could influence 
or could be influenced by glitch events, nor how these might 
rel ate to the quasi-pe riodic spin-down changes observed in 
the lLvne et al.l (|2010D sample. 

Both emission state-changing and glitch activity are 
associated with a wide range of timescales which, for 
each phenomenon, are broadly consistent across discrete 
events in an individual object. For example, the time be- 
tween glitches has a characteristic timescale for each pul- 
sar, and in a few c ases exhibits quasi-periodic behaviour 
IjMelatos et 81.1 12008). The time spent in each emission state 
also shows quasi-periodic fluctuations. Glitch timescales can 
be u sed to probe the physics of the neutron star interior 
(e.g. Ivan Evsden fc Melatosl 120101 ). If glitches and emission 
state switching are be related then perhaps some of these 
timescales are driven by the same physical processes. 

In the MJD 55022 glitch of PSR J0742-2822, the spin 
period of the pulsar increased by Ais = 0.61 piHz, which can 
be considered to be a lower bound on the differential rotation 
frequency between the neutron star crust and its interior. In 
this case the value of Av implies a differential rotation pe- 
riod between the neutron star crust and interior of less than 
~ 19 days. This is about a factor of 5 smaller than the ob- 
served 'typical' periodicity in the state changes, however this 
differential rotation period is valid only prior to the glitch 
and we are only sensitive to state changes in excess of ~ 60 
days. We do not know what the post-glitch differential rota- 
tion period is, but it should be much longer because angular 
momentum has been transferred from the interior to the 
slower-moving crust. Is it possible then that this differential 
rotation is driving the periodicity in the observed i> and pro- 
file shape changes? We note that PSR J0742-2822 has the 
shortest quasi-periodicity of all the state-changing pulsars, 
and exhibits moderate glitches, implying a moderate dif- 
ferential rotation rate. The other pulsars have longer mode 



6 M. J. Keith et al. 

changing timescales, which we could model as a smaller dif- 
ferential rotation rate and therefore may not be expected 
to glitch as rapidly. However, there is currently little direct 
evidence for any link between the differential rotation pe- 
riod and the spin down rate changes. Detailed studies of 
other state-changing pulsars, and high-cadence studies of 
other likely candidates with short period variation would be 
greatly valuable in testing this idea. 
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